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KEY  WORDS:  Beaches;  Coastal  engineering;  Graphs,  Harbors;  Refraction; 
Waves;  WaVej  height 


l!"3 


ABSTRACT:  A  simplified  method  for  plotting  the  refraction  of  waves  as 
they  traverse  water  of  varying  depths  is  presented.  Traditional  methods 
base  the  refraction  diagram  on  a  series  of  bathymetric  contours  spaced 
at  equal  intervals.  This  necessitates  a  relatively  large  number  of 
computations  of  wave  celerity  and  complicated  plotting  methods.  The 
simplified  method  uses  a  variable  contour  interval  so  that  the  ratio 
between  the  celerity  of  the  wave  on  the  shallower  side  of  the  contour 
and  that  on  the  deeper  side  is  a  constant.  Dimensionless  values  for 
the  depth  of  the  contours  related  to  the  wave  period  are  computed  for  the 
celerity  ratio  desired.  The  method  gives  better  resolution  with  fewer 
steps  in  plotting  than  the  traditional  method.  Templates  for  plotting 
the  refraction  according  to  Snell's  Law  are  greatly  simplified  because 
of  the  constant  ratio  of  the  celerities.  Tables  and  a  refraction  template 
are  provided  as  well  as  a  general  discussion  of  wave  refraction.  A 
REFERENCE:  Michel,  John  F.,  "A  Simplified  Method  for  Plotting  Wave 
Refraction,"  Journal  of  the  Waterways,  Harbours  and  Coastal  Engineering 
Division,  ASCE,  Vol. _ ,  No. _ ,  Proc.  Paper  _ 
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A  SIMPLIFIED  METHOD  FOR  PLOTTING  WAVE  REFRACTION* 

By  John  F.  Michel,1  F.  ASCE 

INTRODUCTION 

Gravity  waves  in  water  crossing  bottom  contours  at  an  angle  to  the 
normal  are  bent  or  refracted  toward  the  shallower  contours.  This  is  caused 
by  the  fact  that  the  portion  of  a  wave  in  deeper  water  moves  with  a  greater 
velocity  than  the  portion  in  shallow  water.  Assuming  that  no  energy  is 
transmitted  parallel  to  the  wave  crest,  the  amount  of  energy  between  any  pair 
of  orthogonals  or  lines  drawn  perpendicular  to  the  crest  will  remain  constant. 
Therefore,  the  energy  per  unit  length  of  wave  crest  will  vary  inversely  as  the 
distance  between  the  orthogonais.  Since  the  energy  of  a  unit  length  of  crest 
is  proportional  to  the  square  of  the  wave  height,  the  wave  height  as  modified 
by  refraction  will  vary  inversely  as  the  square  root  of  the  distance  between 
orthogonals.  In  addition  to  the  effect  of  refraction  it  should  be  realized 
that  shoaling  also  causes  modification  of  the  wave  height.  This  will  be 
discussed  in  more  detail  later. 

Refraction  diagrams  such  as  the  one  shown  on  Figure  1  are  extremely 
useful  in  predicting  wave  height  and  angle  of  approach  in  the  nearshore 
region.  The  construction  of  such  a  diagram  is  based  on  the  relationship 
between  the  phase  velocity  or  celerity  C  and  the  depth  h  for  a  wave  having 
a  length  L,  which  is  (2) 

aThe  manuscript  of  this  paper  was  submitted  to  ASCE  on  10  December 
1973  for  publication  in  the  Journal  of  Waterways,  Harbours  and  Coastal 
Engineering  Division. 
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FIGURE  A-6  REFRACTION  DIAGRAM 


A-47 


FIG.  1. — Refraction  Diagram  (Ref.  1,  p.  A-47) 
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C 


tanh 


¥) W2 


(1) 


and  the  relationship  of  length  to  celerity  and  period  T  which  can  be  shown 
tQ  be  constant  for  any  train  of  waves  where 

L  =  CT  or  |  «  Constant  ^ 

Thus  the  ratio  between  the  length  and  celerity  of  a  wave  in  depths  h^  and 
h2  will  be: 


L 

L 


1 

2 


2rrh, 


tanh 


2nh- 


tanh 


(3) 


h  2nh 

For  ratios  of  —  >  .5,  tanh  -7—  is  nearly  unity  so  that  for  deep  water 
L 

conditions  the  relationships  denoted  by  a  subscript  "o"  become: 


C 

o 


fM  1/2 

,  2n  / 


Si 

2n 


(4) 


and 


(5) 


From  this  it  is  seen  that 


C 

C 


tanh 


2rrh 

L 


L 

L. 


(6) 


Tables  (5)  are  available  which  give  values  of  tanh  — for  specified  values 

of  7  and  ^  which  Is  easier  to  compute  since  L  is  independent  of  h. 

L  L  ® 
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It  is  shown  (1)  that  a  wave  crest  tracking  in  a  depth  h,  at  an  angle 
to  a  bathymetric  contour  will  conform  to  Snell's  law  or: 


sin  a„ 


3-  .in  «1 


substituting  Equation  3  in  Equation  7  gives: 

2rrh„ 


tanh 


sin  o. 


tanh  2rrh. 


sin  a. 


(7) 


(8) 


On  the  basis  of  these  relationships,  orthogonals  or  wave  fronts  can 
be  plotted  and  the  wave  heights  related  to  the  deep  water  wave  height  by: 

1/2 

D. 

(9) 


where  H  is  the  wave  height  at  a  given  location,  Hq  is  the  wave  height 

in  deep  water,  bQ  is  the  distance  between  the  orthogonals  in  deep  water, 

b  is  the  distance  between  orthogonals  at  the  given  location  and  D.  or  — , 

an 

o 

is  the  shoaling  coefficient  which  expresses  the  effect  of  depth  on  the 
unrefracted  waves. 

Dj 


■mm 


(10) 


where 


1  + 


4nh 

L 


sinh 


4nh 


(11) 


Values  of  n  or  can  be  computed  or  obtained  from  the  tables (5)  mentioned 


above 


THE  TRADITIONAL  APPROACH 


In  the  traditional  approach  (1)  to  plotting  refraction  diagrams  a 
series  of  bathymetric  contours  with  more  or  less  uniform  spacing  such  as  6 
feet  are  plotted  and  smoothed  and  the  ratios  of  celerities  are  computed. 

Figure  2  typifies  this  approach.  Note  that  for  a  wave  having  a  period  of 
10  seconds  as  specified,  Lq  =  512  feet  so  that  the  table  would  be  carried 
out  to  a  depth  of  255  feet  requiring  43  separate  depths  for  the  6-foot 
interval  shown.  In  actual  practice  as  the  ratio  between  the  celerities  of 
deeper  and  shallower  water  approaches  unity  the  depth  interval  is  increased 
by  leaving  out  some  contours.  The  intervals  however,  are  maintained  in  a 
whole  number  of  feet  or  fathoms. 

Once  the  celerities  at  the  various  depths  have  been  determined,  two 
methods  are  available  for  plotting  the  diagram:  the  wave  front  method  as 
described  by  Weigel  (4)  and  the  more  commonly  used  orthogonal  method  as 
described  by  Johnson  et  al  (3). 

The  Wave  Front  Method 

In  the  wave  front  method,  the  position  of  every  nth  wave  crest  is  plotted 
by  scaling  n  times  the  wavelength  at  the  depth  the  wave  crosses  in  a  direction 
normal  to  the  wave  front.  This  is  done  at  intervals  along  the  crest  until  the 
diagram  is  advanced  to  the  point  of  interest.  Analytic  methods  as  described 
by  Wilson  (6)  utilize  essentially  this  method  for  application  to  digital 
computers . 

The  Orthogonal  Method 

In  this  method  plotting  is  accomplished  by  sketching  a  mid-contour 
between  the  selected  contours  to  serve  as  a  point  of  intersection  in  the 


Part  I  Chapter  I 


TABLE  1-7 


COMPUTATIONS  FOR  VALUES  OF  C|/C2 


T  «  10  seconds 

C,/C2 

C2/Cl 

A 

d 

U 

OP  VCs 

~  Vcd 

(1) 

(feet) 

(2) 

(3) 

(4) 

(5) 

6 

0.01 17 

0.268 

1.40 

0.72 

12 

0.0234 

0.374 

1.21 

0.83 

18 

0.0352 

0.453 

1.14 

0.88 

24 

0.0469 

0.516 

Column  I  gives  depths  corresponding  to  chart  contours.  These  would 
extend  from  6  feet  to  a  depth  equal  to  L<>/2. 

Column  2  is  column  I  divided  by  Lq  corresponding  to  the  determined 
period. 

Column  3:  these  values  may  be  found  In  Table  0-1  of  Appendix  0,  as  a 
function  of  d/Lo*  This  term  is  also  C/Co. 

Column  4  is  the  quotient  of  successive  terms  In  column  3. 

Column  5  is  the  reciprocal  of  column  4. 


FIG.  2 — Computations  for  Refraction  Diagram  (Ref.  1,  Table  1-7,  p.  69) 
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orthogonal.  The  previously  computed  orthogonal  is  extended  to  the  mid¬ 
contour  and  the  angle  cr  of  the  new  orthogonal  is  calculated  using  Snell's 
law  (equation  7).  The  outgoing  orthogonal  is  then  plotted  either  with  a 
protractor  or  a  fairly  complicated  template. 

THE  SIMPLIFIED  METHOD 

In  1969,  the  author,  faced  with  the  problem  of  preparing  some  refraction 

diagrams  in  a  limited  time,  conceived  that  good  resolution  could  be  obtained 

with  less  work  by  maintaining  a  constant  ratio,  R,  between  C  and  C,  so  that 
”  s  d 


R  = 


C 

s 


Cons tant 


(12) 


This  ratio  is  selected  according  to  the  degree  of  resolution  desired.  Values 
of  0.8  or  0.9  yield  good  results  and  are  demonstrated  here.  Once  the  ratio  R 
has  been  selected,  a  universal  table  of  the  depths  of  mid-contours  applicable 
to  any  wave  period  can  be  developed. 

C 

First  values  of  _s  for  each  depth  interval,  are  determined  by: 

C 

o 


f.  -/-w 

Co  (13) 

where  n  is  the  number  of  the  interval.  The  celerity  at  the  mid-contour  is 
then  computed  as  the  average  of  and  Cg  or 


so  that 


mid 


C  +  C . 
s  d 


Cmid  -  \  (1+R)  R(n_1)  -  K 

A  ^ 


(14) 


(15) 


which  is  a  till  a  constant  for  any  value  of  n 


At  this  point,  it  is  well  to  note  that  the  depth  of  the  mid-contour 
is  selected  as  the  depth  corresponding  to  the  average  celerity  over  the 
depth  interval.  The  traditional  approach  on  the  other  hand  utilizes  as 
a  mid-depth  a  contour  placed  by  eye  between  the  computed  depths. 

Since : 


K  = 


Cmid 


2nh 


tanh 


mid 


mid 


mid 

L 


where  L  ^  is  the  wavelength  at  the  mid-depth, 


2nh  . 

mid  ^  - 1  „ 

— -  =  tanh  K 

mid 


(16) 


(17) 


and 


Since 


h 

mid 


i 


iii 

2tt 


(18) 


mid  = 


gT 


K  tanh 

(2n)2 


(19) 


which  is  dimensionless  and  dependent  only  in  the  constant  ratio  and 

hmid  d 

the  wave  period.  A  more  convenient  ratio,  — =—  in  meters  or  feet  can  be 

TZ 

obtained  by  multiplying  the  expression  by  g  in  the  units  desired. 

At  this  stage  values  of  the  shoaling  coefficient  D,  or  — ,  can  be  computed  for 

Q  H 

h  ° 

each  selected  ratio  of  —  • 

Tables  1  and  2  illustrate  the  relationship  described  above  for  two 
C 


values  of 


which  will  cover  a  reasonable  range  of  wave  periods. 


TABLE  2:  COMPUTATIONS  FOR  MID-CONTOURS  WITH  C  /C,  *  0.90 

s  d 


o 

o 

0} 

o 

(1) 

o 

o 

*H 

B 

O 

(2) 

CM 

H 

00 

•o 

•w 

(3) 

CM-s 

PJ 

O 

u 

(D 

CM  0) 

H  \ 

w 

•O 

•H  Q) 

B  u 

M  QJ 

B 

(4r 

hmld/T 

2 

(feet/second  ) 

£  H/H' 

0 

1.0000 

.9500 

■ 

.4320 

1.4173 

.94 

.9000 

.8550 

HgSl 

.2705 

.8874 

.91 

.8100 

.7695 

.1947 

.6387 

.92 

.7290 

.6926 

HI 

.1466 

.4810 

.94 

.6561 

.6233 

.0115 

.1130 

.3708 

.97 

.5905 

.5610 

.0090 

.0883 

.2898 

1.00 

.5314 

.5049 

.0071 

.0697 

.2285 

1.04 

.4783 

.4544 

.0056 

.0553 

.1814 

1.09 

.4305 

.4089 

.0045 

.0441 

.1447 

1.14 

.  3S74 

.3680 

.0036 

.0353 

.1158 

1.19 

.3487 

.3312 

.0029 

.0283 

.0929 

1.25 

.3138 

.2981 

.0023 

.0228 

.0746 

1.32 

.2824 

.2683 

.0019 

.0183 

.0601 

1.38 

.2542 

.2415 

.0015 

.0148 

.0484 

1.45 

.2288 

.2173 

.0012 

.0119 

.0391 

1.53 

.2059 

.1956 

.0010 

.0096 

.0316 

1.61 
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Actually  only  the  third  column  is  necessary  for  the  computation  of  the 
depth  of  the  mid-contour,  the  fourth  and  fifth  columns  having  been  added 
for  convenience.  The  first  two  columns  are  shown  only  to  demonstrate  the 
rationale  of  the  method.  The  sixth  column  provides  ready  reference  to 
the  shoalling  coefficient.  These  tables  are  universal  in  application  and 
are  sufficiently  brief  to  be  placed  on  a  refraction  template  as  will  be 
seen  later. 

In  utilizing  Table  1  for  instance,  for  a  wave  having  a  period  of  10 
seconds,  the  depths  of  the  mid-contours  in  feet  could  be:  108,  53,  31,  19, 
12,  7  and  5.  These  depths  would  be  plotted  on  the  chart  to  be  used.  If 
Table  2  were  used  for  better  resolution,  the  depths  would  be  142,  89,  64, 
48,  37,  29,  23,  18,  14,  12,  9,  7,  6  and  5.  These  14  depths  can  be  compared 
to  Figure  2  which  would  require  43  depths  to  be  plotted. 

When  the  mid-contours  have  been  plotted,  either  the  wave  front  method 
or  the  orthogonal  method  can  be  used  with  great  facility  to  plot  the  re¬ 
fraction  diagram. 

The  Wave  Front  Method 

In  using  the  wave  front  method,  a  set  of  proportional  dividers  set  to 
C 

the  ratio  —  can  be  used  to  plot  the  positions  of  the  nth  successive  wave 
Cd 

fronts  without  the  necessity  for  any  measuring  or  scaling  except  for  the 
crests  in  deep  water  in  case  it  is  desired  to  show  actual  crest  position 
instead  of  simply  the  shape  of  the  wave  front. 

Although  it  has  not  been  tried  this  method  would  appear  to  be  well 
adapted  to  numerical  analysis,  by  digital  computer,  and  its  simplicity 
should  greatly  reduce  computation  time. 
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The  Orthogonal  Method 

In  utilizing  the  otthogonal  method  the  constant  relationship  between 
sin  and  sin  permits  construction  of  a  very  simple  template.  Table 

c 

g 

3  shows  the  relationships  between  orthogonal  angles,  for  ratios  of  r-  ■  0.9 

Cd 

and  0.8. 

C 

g 

An  example  of  such  a  template  for  —  =  0.9  is  shown  in  Figure  3. 

a 

Since  the  template  has  graduations  in  both  sides,  there  is  no  need  to  have 
it  on  transparent  material. 

In  using  the  template,  the  depths  of  the  mid-contour  are  determined 
from  Table  2  or  the  slightly  abbreviated  table  on  the  template.  When  the 
mid-contours  have  been  plotted  on  the  chart  the  template  is  used  to  plot  the 
orthogonals  in  the  following  manner  as  demonstrated  by  Figure  4: 

1.  The  template  is  placed  so  that  the  tangent  line  is  on  the  tangent 
to  the  mid-contour  at  the  point  where  the  incoming  orthogonal  intersects  it. 
This  places  the  baseline  of  the  template  on  the  normal  to  the  contour  which 
is  drawn  using  the  base  as  a  straight  edge. 

2.  The  angle  of  the  incoming  orthogonal  with  the  normal  is  noted  on  the 
scale  for  deeper  or  shallower  water  as  the  case  may  be. 

3.  The  template  is  then  rotated  so  that  the  angle  on  the  other  side  of 
the  template  corresponding  to  the  incoming  angle  is  on  the  normal  line.  The 
new  orthogonal  is  then  drawn  using  the  base  of  the  template  as  a  straight 
edge. 


TABLE  3:  SIMPLIFIED  METHOD  COMPUTATION  FOR  REFRACTION  TEMPLATES 


a. 

<Vcd  ■  •» 

a 

C  /C,  -  .8 
s  a 

a 

d 

s 

3 

(1) 

(2) 

(3) 

degrees 

degrees 

degrees 

0 

0 

0 

10 

9.0 

8.0 

20 

17.9 

15.9 

30 

26.7 

23.6 

40 

35.3 

30.9 

50 

43.6 

37.8 

60 

51.2 

43.9 

70 

57.7 

48.7 

80 

62.4 

52.0 

90 

64.2 

53.1 

2.  Place  template  so  that  normal  line  is  on  same 

GRADUATION  ON  ARC  ON  THE  OTHER  SIDE  AS  NOTED  ABOVE  FOR  INCOMING 
ORTHOGONAL.  DRAW  OUTGOING  ORTHOGONAL^ ALONG  BASELINE. 


Cs/Cd=0.  90 


SHALLOW  SIDE  h 

DEEP  SIDE 

Z 

III 

bJ 

z 
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Z 

J 

< 

h 

BASE 

LINE 

FIG.  3 — Refraction  Teaplate 
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SUMMARY 

The  method  of  plotting  wave  refraction  diagrams  shown  here  is  extremely 
efficient  and  simple  to  use  because  the  ratio  of  shallower  to  deeper  water 
celerity  is  maintained  constant.  The  use  of  universal  tables  and  simple 
templates  replaces  a  considerable  amount  of  computation  and  tedious 
plotting  utilized  in  the  traditional  method.  It  also  provides  better 
resolution  by  utilizing  mid-contours  that  are  closer  together  in  shallower 
water.  The  simplified  method  can  be  used  either  with  the  wave  front  method 
or  the  orthogonal  method  of  plotting.  It  may  possibly  facilitate  in  the 
same  manner  the  utilization  of  numerical  methods  with  digital  computers. 


REFERENCES 


1.  Corps  of  Engineers,  U  .S .  Army,  Coastal  Engineering  Research  Center 
"Shore  Protection,  Planning  and  Design",  Technical  Report  No.  4,  June 
1966,  pp.  63-72. 

2.  Ippen ,  A.  T.,  "Estuary  and  Coastline  Hydrodynamics",  McGraw-Hill,  1966, 
pp.  23-26. 

3.  Johnson,  J.  W.,  0  Brien,  M.P.,  and  Issacs ,  J.D.,  "Graphical  Construction 
of  Wave  Refraction  Diagrams",  H.O.  Pub.  No.  60S,  U.S.N.  Hydrographic 
Office,  January,  1946. 

4.  Wiegel,  R.  L.,  "Oceanographical  Engineering",  Prentice  Hall,  1964. 
pp.  163-168. 

5.  Wiegel,  R.  L.,  "Gravity  Waves,  Tables  of  Functions",  Coastal  Engineering 
Research  Council,  A.S.C.E.,  Reprinted  1967. 

6.  Wilson,  W.  S.,  "A  Method  for  Calculating  and  Plotting  Surface  Wave  Rays", 
Technical  Memo  No.  17,  U.S.  Army  Coastal  Engineering  Research  Center, 
February  1966. 


* 


18 


APPENDIX  -  NOTATION 


The  following  symbols  are  used  in  this  paper: 


b 

b 

o 

C 

C1’C2 

VCd 

Cmid 

C 

o 

Dd 

g 

h 

Vh2 

*\nid 

H 

H  ' 
o 

H/H  ' 
o 

K 

Ll-L2 

L 

n 

R 

T 

Of 


distance  between  orthogonals; 

distance  between  orthogonals  in  deep  water; 

celerity,  a  phase  velocity  of  a  wave; 
celerity  of  a  wave  in  depths  h^  and  h^; 

celerity  of  a  wave  in  shallower  or  deeper  water; 

celerity  of  a  wave  at  mid-contour; 

celerity  of  a  wave  in  deep  water; 

shoaling  coefficient  or  ratio  of  unrefracted  wave  height 

to  deep  water  wave  height; 

gravitational  constant; 

depth  of  water  to  still  water  level; 

depths  corresponding  to  C^.C^; 

the  depth  at  a  mid-contour; 

wave  height  from  trough  to  crest; 
unrefracted  wave  height  in  deep  water; 

see  above; 

a  constant  ratio  of  C  .  ,/C  ; 

mid  o 

wave  length  in  depths  h^.h^; 

wave  length  measured  from  crest  to  crest; 

the  number  of  a  depth  interval; 

a  constant  ratio  of  C  / C,; 

s  d 

the  period  of  a  wave; 

the  angle  between  a  wave  crest  and  a  bathymetric  contour; 
angles  between  wave  crests  and  contours  in  shallower  or 
deeper  water. 


APPENDIX  -  COMPUTER  PROGRAM 
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